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FFLO or Majorana superfluids: Quantum phases of fermionic cold atoms in spin-orbit 

coupled optical lattices 

Chunlei Qu 1 , Ming Gong 1 , and Chuanwei Zhang^ 
1 Department of Physics, The University of Texas at Dallas, Richardson, TX, 75080 USA 

The recent experimental realization of spin-orbit coupling (SOC) for ultra-cold atoms opens a 
completely new avenue for exploring new quantum matter. In experiments, the SOC is implemented 
simultaneously with a Zeeman field. Such spin-orbit coupled Fermi gases are predicted to support 
Majorana fermions with non-Abelian exchange statistics in one dimension (ID). However, as shown 
in recent theory and experiments for ID spin-imbalanced Fermi gases, the Zeeman field can lead 
to the long-sought Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superfluids with non-zero momentum 
Cooper pairings, in contrast to the zero momentum pairing in Majorana superfluids. Therefore a 
natural question to ask is which phase, FFLO or Majorana superfluids, will survive in spin-orbit 
coupled Fermi gases in the presence of a large Zeeman field. In this paper, we address this question 
by studying the mean field quantum phases of ID spin-orbit coupled fermionic cold atom optical 
lattices. 

PACS numbers: 03.75.Ss, 67.85.-d, 74.20.Fg 



I. INTRODUCTION 

Spin-orbit (SO) coupling plays an important role in 
many important condensed matter phenomena p], [2j, 
ranging from spintronics to topological insulators. The 
recent experimental breakthrough in realizing SO cou- 
pling in ultra-cold Bose and Fermi gases [3-7] provides 
a new platform for engineering many new many-body 
quantum matters [8|- In the experiments, the SO cou- 
pling is realized together with a Zeeman field. It is well 
known that such SO coupling and Zeeman field, together 
with the 5-wave superfluid pairing in degenerate Fermi 
gases, can support zero-energy Majorana fermions [9j 
with non-Abelian exchange statistics when the Zeeman 
field is beyond a certain critical value [10|, [ll| . In solid 
state, the same ingredient has been realized using a het- 
erostructure composed of a semiconductor nanowire (or 
thin film) with strong SO coupling, an s-wave supercon- 
ductor, and a magnetic field (or a magnetic insulator) 
|12-[l9j. Important experimental progresses have been 
made along this direction |20l-[23j , where some signatures 
which may be related with Majorana fermions have been 
observed. In the solid state heterostructure, the s-wave 
pairing is induced to the semiconductor through proxim- 
ity effects [H IH . 

In degenerate fermi gases, however, as observed in 
both theory and experiments [26|-[28j, the presence of a 
large Zeeman field (realized by the spin population im- 
balance) can induce non-zero momentum Cooper pairings 
between atoms, i.e., FFLO phases J29J-|31|, especially in 
low-dimensional Fermi gases. Such FFLO phases may 
not support Majorana fermions. Therefore it is natu- 
ral to ask whether the FFLO superfluids |32l-[39| with 
non-zero momentum Cooper pairs or the Majorana su- 
perfluids with zero momentum Cooper pairs will survive 
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in the presence of SO coupling and a large Zeeman field. 
This question becomes especially important because of 
the recent experimental realization of SO coupled Fermi 
gases ja, LZ|, which makes the observation of Majorana 
fermions in cold atomic systems tantalizingly close. The 
cold atom system may be a better platform for the ob- 
servation of Majorana fermions because of the lack of 
disorder and impurity [36|, |40| - |45J , an issue that has led 
to intensive debate in the condensed matter community 
on the zero-bias peak signature of Majorana fermions in 
recent transport experiments [20|, Um • 

In this paper, we address the competition between 
FFLO and Majorana superfluids by studying the quan- 
tum phases of spin-imbalanced Fermi gases in spin-orbit 
coupled optical lattices. Because of the fact that the 
experimentally realized SO coupling is one-dimensional 
(ID) and the natural dimension requirement for the real- 
ization of Majorana fermions in this system, we consider 
ID SO coupled optical lattices (similar to ID nanowires) 
and investigate the quantum phase diagram at zero tem- 
perature using the mean field theory. The quantum 
phases of Fermi atoms are obtained by self-consistently 
solving the corresponding Bogoliubov-de Gennes (BdG) 
equation. Without SO coupling, there are no Majorana 
superfluids, and FFLO superfluids appear in the large 
Zeeman field region. The SO coupling enhances the Ma- 
jorana superfluid phase while suppresses the FFLO su- 
perfluid phase. Majorana and FFLO superfluids exist for 
different filling factors (i.e., in different chemical poten- 
tial region). We characterize different quantum phases 
by visualizing their real space superfluid order param- 
eters, density distributions, and Majorana zero energy 
wavefunctions. The effects of the harmonic trap are also 
discussed. 

The rest of the paper is organized as follows. In Sec. 
HH we present the BdG equation for describing Fermi 
atoms in spin-orbit coupled optical lattices. The sym- 
metry of the BdG equation and its consequence to the 
phase diagram is discussed. In Sec. IIII| we present the 



main numerical results obtained by self-consistently solv- 
ing the BdG equation. We discuss the phase diagram, the 
characterization of various phases, and the effects of the 
harmonic trap and Hartree shift. Sec. IIVI consists of 
discussion and conclusion. 



where the chemical potential fli a = ji + U(hia) becomes 
site dependent, a = —a. The superfluid pair potential 
is defined as A^ = — [Z(c^c^). Using the Bogoliubov 
transformation Ci a = J2 n (u™ a j n — av^jl) , we obtain the 
BdG equation 



II. HAMILTONIAN AND SYMMETRY 

We consider a ID degenerate Fermi gas with Zeeman 
field and SO coupling in an optical lattice. The dynamics 
of this system can be described by the standard tight- 
binding Hamiltonian 



H = U + Hz + Hso, 



(1) 



where the first term is the usual spin- 1/2 Fermi-Hubbard 
model in an optical lattice, 

Uo = -t y^(4 a Q+ia + H.c.) - /i^hia -U^2 ft *t™*+> 

ia ia i 

(2) 
t is the hopping amplitude, \i is the chemical potential, 
and U is the contact interaction. The second and third 
terms are the Zeeman field 

i 

and SO coupling 

i 

Such type of SO coupling and Zeeman field have been re- 
alized in recent experiments for both bosons and fermions 
[3|-l7J]. In experiments, the hopping amplitude, Zeeman 
field, SO coupling strength, and contact interactions may 
be tuned independently. In our numerical simulation we 
set t — 1 throughout this work. All other energies are 
scaled by t. The total length of the ID optical lattice is 
chosen as N = 100, which is long enough to ensure that 
the coupling between the two ends is vanishing small. 
An open boundary condition is used to obtain the zero- 
energy Majorana fermions at the two ends of the ID op- 
tical lattices in the topological superfluid regime. 

As the first approach for understanding the quantum 
phases of such ID SO coupled optical lattices, we con- 
sider the standard mean field theory. We decouple the 
interaction term in Hq using 



-Un^fiii = Aic]+c. 
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Notice that here we have taken into account the Hartree 
shift term, which has quantitative effects on our results. 
The effective Hamiltonian reads as 

H eS =~tJ2 Z)(4*+l* + H.C.) - J2 hAaCia 

+ J2( A 4Al + A *<M*t) + n z + Hso, (3) 
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with the Fermi-Dirac distribution f(E) = 1/ (l + e E/T ). 
The BdG equation (j4]) should be solved self-consistently 
with the order parameter equation ([5]) and the particle 
number equation (J6j) for the average number of atoms per 
lattice site n = J2 { a (hi a ) /N. Here we denote n as the 
filling factor for convenience (note that generally the fill- 
ing factor is defined as v = n/2 in the literature). In our 
simulation, we take the attractive interaction strength 
U = 4.5 and the temperature T = 0. 

The original Hamiltonian Eq. (pQ) contains no imagi- 
nary part, therefore the wavefunction in Eq. (jlj) can be 
made real and the order parameter in Eq. ([3]) is also 
real. The topological symmetry of the effective model 
belongs to BDI class, with characteristic index Z, in- 
stead of the BdG D-class with Z 2 index [13, E3|. The 
mean field Hamiltonian ([3]) still preserves the basic sym- 
metry properties of Eq. (pQ). Consider the particle- 



hole operation C 



dt\ c -i = (-i)MT 



we have 



CH"(/i)C _1 = H(—/jl) when the Hartree shift term is ig- 
nored and CH"(/i)C _1 = H(—U — n) when the Hartree 
shift term is included. Thus the spectrum should be sym- 
metric about \± — and /i = — ^ respectively, which is 
confirmed in our numerical results (see Fig. 1 and Fig. 
6). However, in the presence of a trapping potential, 
the chemical potential becomes site dependent, thus the 
particle-hole symmetry C is broken and the band struc- 
ture is no longer symmetric about ji = (see Fig. 5). 
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FIG. 1. (Color online) Phase diagram of ID SO coupled opti- 
cal lattices as a function of Zeeman field h and filling factor n 
(a,c) or chemical potential \i (b,d). First row: the spin orbit 
coupling a = 0.5, second row: a = 1.0. Five different phases 
are identified in each phase diagram: normal BCS superfluid 
(NS), topological superfluid(TS), FFLO, normal gas (NG), 
and insulator phase (Ins). 
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FIG. 2. (Color online) Representative density of states for (a) 
insulator phase, (b) normal BCS superfluid phase, (c) topo- 
logical superfluid phase, and (d) FFLO phase. The insets 
show the order parameters in each phase. The corresponding 
phase points are marked by the plus signs in Fig. [TJd. Note 
that there is a zero-energy peak in the DOS in the topological 
superfluid phase as shown in (c). 



III. 



PHASE DIAGRAM: MAJORANA VERSUS 
FFLO SUPERFLUIDS 



We self-consistently solve the BdG equations (j4|5|6p 
with an open boundary condition to obtain the phase 
diagram. The SO coupled optical lattice supports sev- 
eral different phases: normal BCS superfluid (NS) with 
A / and all non-zero eigenstates; topological super- 
fluids with A / and zero-energy Major ana fermions 
located at two ends of the lattice; FFLO phase with os- 
cillating A and magnetization; insulator phase (Ins) with 
integer filling factor and finite energy gaps; and normal 
gas (NG) phase without pairing and energy gap. We first 
study the phases without Hartree shift, and address the 
role of Hartree shift at the end of the section. 



A. Phase diagram without Hartree shift 

Our numerical results are presented in Fig. [1] for two 
different sets of SO coupling strength. In Figs. [1] (a) and 
(c) we plot the phase diagram in the h — n plane and in 
Figs. IU(b) and (d) we plot the results in the h — fi plane. 
We see the phase diagram is symmetric around n = 1 or 
\i = 0, as discussed in the previous section. When the 
Zeeman field h is very small, the system favors the nor- 
mal BCS superfluids. With increasing Zeeman field /z, 
topological superfluids and FFLO phases emerge as the 
ground states of the system for different filling factors. 
The FFLO phase is more likely to be observed around 



the integer filling factor n — 1. When the Zeeman field 
becomes even larger, where only atoms with one type of 
spin can stay in each lattice site, the insulator phase de- 
velops with the filling factor n = 1. The insulator phase 
can also be found when \i is too large (fully occupied 
band) or to small (empty band). The topological phase 
and associated Major ana fermions emerge with fractional 
filling factors n. By comparing the phase region for dif- 
ferent SO coupling, we find that the SO coupling en- 
hances the topological superfluid phase and suppresses 
the FFLO phase. Note that without SO coupling, there 
is only FFLO phase, and no topological superfluid phase 
[H,[2Z|. 

Different quantum phases in Fig. [1] can be char- 
acterized by different density of states (DOS) p(E) 



Y:[\u ia \ 2 S(E-E n )- 



^S(E + E n )] and superfluid or- 



der parameter A (x), as shown in Fig. O In the insulator 
phase (Fig. [2^), the order parameter A = and there is 
an energy gap for excitations. In the normal BCS super- 
fluid phase (Fig. [2}d), the order parameter is non-zero and 
there is a superfluid gap around E = 0. In the topological 
superfluid phase (Fig. [2t), a zero energy peak appears in 
the DOS which corresponds to the zero energy Majorana 
state. Note that for the normal BCS superfluid phase the 
order parameter has a strong oscillation near each end, 
while for the topological superfluid, the order parame- 
ters is a monotonic function of the sites. Similar features 
have always been found for these two different phases for 
different parameters. The FFLO phase (Fig. [2]i) is char- 
acterized by the spatially oscillating order parameter and 
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FIG. 3. (Color online) (a) Plot of the order parameter and the 
lowest two eigenenergies E\ and E2 as a function of Zeeman 
field for the transition from normal BCS superfluids to topo- 
logical superfluids. a = 1.0, fi = —2.0. (b,c) The Majorana 
zero energy state wavefunction. 



4 


Ins 


2 


f5c::is::7f: NS 

NG \ / " 


W 


- r^IIX^CII^FFLO 


-2 


- NG \ / " 




-4 


" (a) . ,1ns ....." 




FIG. 4. (Color online) Single particle band structure of the 
spin-orbit coupled optical lattices, a = 1.0, h = 1.7 (corre- 
sponds to the dashed line in Fig. [TJl). Corresponding chem- 
ical potential regions for different phases are identified, (b) 
Plot of /x as a function of Zeeman field for fixed filling factor 
n — 0.25, 0.5, 0.75, 1.0. The regions between the two circle 
symbols on each line label the topological superfluid phase 
(for n = 0.25, 0.5) or the FFLO phase (for n = 0.75, 1.0). 



the excitations are always gapless. 

The emergence of Majorana zero energy state at the 
ends of the SO coupled optical lattices can be clearly seen 
in Fig. [3j In Fig. [3] we plot the bulk order parameter, 
the first and second non- negative eigenvalues, denoted 
by Ei and E2 , of the BdG equation as a function of Zee- 
man field. We see that in the normal BCS superfluid 
regime E\ is always equal to £"2, and smaller than the 
order parameter. The energy gap not equal to the or- 
der parameter is a unique feature for spin-orbit coupled 



systems. When h approach 1.0 (h > A), we observe 
a sudden jump of the order parameter and the system 
enters the topological superfluid regime, where E\ = 0, 
and E^ increases and reaches the maximum value 0.5 at 
h ~ 1.1. When h further increases, E2 ~ A gradually 
decreases and becomes zero when the system enters the 
normal gas phase. In the topological superfluid phase, E2 
is the minimum energy gap that protects the topological 
zero energy Majorana state. In Fig. [3]o and c, we plot 
the wavefunction of the zero energy state, from which we 
confirm that the corresponding Bogoliubov quasiparticle 
operators 7 n = J2ia( u ia c iv + v ia c L) at tne two ends for 
the zero energy states satisfy 7q =70, that is, Majorana 
operators. 

Different phases in different parameter regions can 
be intuitively understood from the single particle band 
structure, as shown in Fig. \$k for a = 1, h = 1.7 (cor- 
responds to the dashed line in Fig. [T]i). The chemical 
potential regions for different phases are identified in the 
figure. Comparing Fig. [T]i and Fig. SH, we see the 
topological superfluid phase appears when the chemical 
potential cuts only a single band, while the FFLO phase 
appears mainly around /i = (filling factor n = 1) and 
cuts two bands. The insulator phase appears when two 
bands are either both fully occupied or empty. Normal 
superfluid or normal gas phases also appear when a single 
band is occupied with either small or large filling factors. 

Because of the SO coupling, the spin polarization is 
not a conserved quantity anymore, which is different from 
the spin-imbalanced Fermi gases in the literature. How- 
ever, the filling factor can still be controlled precisely in 
experiments. With increasing Zeeman field, the chem- 
ical potential changes for a fixed filling factor, leading 
to the transition between different phases. In Fig. HJd 
we plot the chemical potential as a function of Zeeman 
field for different filling factor n. Generally when n < 0.7 
or n > 1.3, we find \i changes monotonically as a func- 
tion of the Zeeman field. However, in the case n = 0.75, 
the chemical potential does not change when h < 1 and 
then increases slightly when the system enters the FFLO 
phase regime. For n — 1, the chemical potential is in- 
dependent of the filling factor due to the particle-hole 
symmetry. Notice that n is not a unique function of the 
chemical potential, therefore in Fig. [1] we see that the 
FFLO phase is very large in the h — n plane but becomes 
much smaller in the h — \i plane because the mapping is 
a very complex function. 



B. Phase diagram in a harmonic trap 

In a realistic experiment, a harmonic trapping poten- 
tial V(x) = ^uj^(i — L c ) 2 exists, where u x is the trapping 
frequency and L c is the center of the lattice. The effects 
of the trapping potential on the filling factor and order 
parameter are shown in Fig. [5fa) and (b). With increas- 
ing trapping frequency, the ultracold atoms are forced to 
the center of the trap which has a lower potential. With 




FIG. 5. Phase diagram in the presence of a harmonic trap. 
The order parameter profiles (a) and atom density distribu- 
tions (b) for different trapping frequencies: uj x = 0.0 (Red 
dotted line), 0.05 (Orange dashed line), 0.15 (Green solid 
line), and 0.2 (Blue dash-dotted line), (c) and (d) are phase 
diagram as a function of Zeeman field h and filling factor n 
or chemical potential ji. 



a high trapping frequency, a shell structure is developed, 
where superfluid, insulator, and normal gas phases ap- 
pear in different regions of the harmonic trap. In the 
superfluid regime in the trap, topological superfluid or 
FFLO phase may develop for different parameters. In 
Fig. [5] (c) and (d) we plot the phase diagram for a fixed 
trapping potential. Notice that there is always a mixture 
of different phases in the trap, therefore we only iden- 
tify three different cases in our plot: normal gas for the 
whole lattice, topological superfluids (with zero-energy 
Majorana fermions) and other superfluids (normal BCS 
superfluids or FFLO) in certain part of the lattice. In 
the topological superfluids, the Majorana wave function 
does not localize at the two ends of the trap, but in cer- 
tain middle region of the trapping potential [37J. The 
position of the Majorana fermion changes with the trap- 
ping frequency. With the trapping potential, the effec- 
tive chemical potential is site-dependent, therefore the 
particle-hole symmetry is broken and the phase diagram 
is not symmetric with respect to n — 1 (Fig. [5] (c)) or 
M = (Fig. 0(d)). 



C. Effect of Hartree shift 

We also study the effect of the Hartree shift on the 
phase diagram. When the Hartree shift is included, the 
local chemical potential and Zeeman field are both mod- 
ified as seen from jl^ = /i + U{riio). If we define (rti) 
and (rrii) as the local particle number and magnetization 
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FIG. 6. Phase diagram with the inclusion of the Hartree shift. 
a — 0.5. All other parameters and notations are the same as 
that in Fig. [Q 



mi — ni\ — n%\,, then the modified chemical potential and 
Zeeman field with including the Hartree shift term will 
be Hi = ji + U (rii) and hi = h — |(mi). The numerical 
results are presented in Fig. [6] for a direct comparison 
with the results without Hartree shift (Fig. [1] (c) and 
(d)). We see the phase diagram is still qualitatively the 
same except that now the phase diagram is symmetric for 
\i = —U/2 as discussed in the above. Since the Zeeman 
field hi depends strongly on the local magnetization rrii 
and the FFLO phase is a competition between superfluid 
and magnetization, it is hard to numerically determine 
FFLO phases near the phase boundary between FFLO 
and topological superfluid phases. This is because the 
total free energy becomes extremely complex as a func- 
tion of order parameters and other quantities, and most 
solutions we find correspond to excited states, instead of 
global minimum of the free energy. Thus the boundary 
between topological superfluid and FFLO phase cannot 
be determined precisely. 



IV. DISCUSSION AND CONCLUSION 

So far we have considered the quantum phases of spin- 
orbit coupled Fermi gases in ID optical lattices using the 
mean field theory. In a truly ID system, the quantum 
fluctuation needs to be taken into account and further 
investigation of quantum phases is needed using exact 
method such as DMRG. The quantum fluctuation may 
also be suppressed by considering a 3D optical lattice 
with weak tunneling along the y and z direction. Such 
tunneling may weaken the quantum fluctuation, but does 
not change the mean field phase diagram qualitatively, as 
confirmed in our numerical calculation. Another method 
for suppressing quantum fluctuation is to consider multi- 
bands along the y and z directions [47j. Nevertheless, our 
results provides an intuitive and first step understanding 
of the competition between FFLO and Majorana super- 
fluids in the ID SO coupled lattice system, and may also 
have application for the solid state nanowire heterostruc- 
ture where similar physics exists. 



In summary, in this paper, we address the question 
that which phases, FFLO or Majorana superfTuids, will 
survive in spin-orbit coupled Fermi gases in the presence 
of a large Zeeman field through studying the mean field 
quantum phases of ID spin-orbit coupled optical lattices. 
In the optical lattice system each site can host at most 
2 fermions, making the system host plenty of phases de- 
pending on the filling factor and the Zeeman field, which 



are quite different from the free space results. At a finite 
Zeeman field we observe the strong competition between 
topological superfluid phase and FFLO phase. The SO 
coupling enhances the topological superfluid phase while 
suppresses the FFLO phase. 
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